Why is Understanding Glassy Polymer Mechanics So Difficult? 
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In this Perspective, I describe recent work on systems in which the traditional distinctions between 
(i) unentangled vs. well-entangled systems and (ii) melts vs. glasses seem least useful, and argue 
for the broader use in glassy polymer mechanics of two more dichotomies: systems which possess 
(iii) unary vs. binary and (iv) cooperative vs. noncooperative relaxation dynamics. I discuss the 
applicability of (iii-iv) to understanding the functional form of strain hardening. Results from 
molecular dynamics simulations show that the "dramatic" hardening observed in densely entangled 
systems is associated with a crossover from unary, noncooperative to binary, cooperative relaxation 
as strain increases; chains stretch between entanglement points, altering the character of local 
plasticity. Promising approaches for future research along these lines are discussed. 




Strain 

FIG. 1: Schematic of stress-strain curves for ductile polymer 
glasses. Dramatic hardening coincides with the increase in 
slope at large strains. In brittle systems, fracture intervenes 
at lower strains because strain hardening is insufficient to sta- 
bilize the material against post-yield strain localization 1]. 

One of the reasons dense polymeric systems so interest 
physical scientists, apart from their ubiquity and utility, 
is the wide range of energy, length, and time (EL&T) 
scales controlling their properties. Figure [TJ a schematic 
depiction of typical stress-strain curves for ductile poly- 
mer glasses, illustrates the mechanical consequences of 
this range. Undeformed systems occupy low-lying regions 
of a rugged free energy landscape. In the linear elastic 
regime, systems remain near initial free energy minima, 
and stress is controlled by local forces at the Kuhn scale 
or below. Yield occurs when energetic barriers to seg- 
mental rearrangements are overcome; the resulting in- 
crease in local mobility produces strain softening. In the 
plastic flow regime, the stress a = dW/de is relatively 
constant. Strain hardening begins when dW/de must in- 
crease to drive further segmental rearrangements while 
maintaining chain connectivity. This increase becomes 
more dramatic as the scale over which chains are oriented 
approaches that of the entanglement mesh. Finally, frac- 
ture occurs when cohesive forces, either primary covalent 
bonds or secondary nonbonded interactions, no longer 
suffice to maintain material integrity. 

The industrial importance of understanding strength 



and failure of polymeric materials has made quantita- 
tively predicting the entire range of mechanical response 
shown in Fig. [T] one of the main goals of physical polymer 
science. While great progress has been made in recent 
years towards understanding phenomena at strains up to 
and including the early stages of strain hardening within 
a single framework, a coherent theoretical picture includ- 
ing dramatic hardening and fracture remains elusive. In 
this Perspective I discuss some of the reasons why this is 
so, and outline possible strategies for moving forward. 

Polymer dynamics and mechanics are often cast in 
terms of polar dichotomies; a given system is classified as 
belonging to one "pole" or the other. The two most com- 
monly employed dichotomies for bulk amorphous systems 
(at least those relevant to Fig. [T]) are those between: (i) 
unentangled vs. well-entangled systems, and (ii) melts 
far above Tg vs. glasses far below Tg. Dynamical behavior 
in polar limits (e. g. unentangled melts 0) is amenable 
to relatively simple theoretical treatment. Because the 
EL&T scales controlling relaxation are well separated, 
one may focus on a dominant scale and then treat relax- 
ation on that scale. In the context of Fig. [U only glassy 
systems posess a finite plastic flow stress, and only en- 
tangled systems strain harden. Flow and hardening have 
therefore traditionally been explained in terms of local 
plasticity [s'l and the work required to deform the entan- 
glement network [4], respectively. 

Dichotomies (i-ii) have obviously proven of great util- 
ity for understanding polymer mechanics, but for actively 
deforming systems, they are less useful. A classic exam- 
ple is the well known tendency of systems above their qui- 
escent Tg to exhibit glass-like mechanical features when 
deformed sufficiently rapidly. Significant strain harden- 
ing occurs when the product of strain rate e and chain- 
scale relaxation time r exceeds unity. More recent exper- 
imental, theoretical and simulation work [^-0] has fur- 
ther blurred the distinction between melts and glasses 
by showing that systems below their quiescent Tg become 
meltlike on short length scales under active deformation. 
Finite strain rates can reduce the segmental relaxation 
time Ta by orders of magnitude; tTa drops well below 
unity at yield. This drop is partially reversed in the 
strain hardening regime, and in all cases, Tq increases 
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dramatically when deformation is ceased. Clearly the 
terms "meltlike" and "glassy" poorly characterize such 
nonlinear behavior. Categorizing systems as either "en- 
tangled" or "unentangled" can be similarly complicated, 
because entanglements (unlike crosslinks) may be short- 
lived or long-lived compared to experimental time scales, 
and their relaxation is also altered by active deformation. 
These developments have helped clarify why developing 
a coherent theoretical framework for predicting the entire 
range of behavior depicted in Fig. [T] is so difficult. 

How, then, to proceed? The great advances made 
by conceptualizing mechanical properties in terms of di- 
chotomies (i-ii) suggest that looking for additional ones 
is a useful strategy. It seems to me that two potentially 
very useful dichotomies for improving our understanding 
of glassy polymer mechanics are: (iii) unary vs. binary 
relaxation and (iv) noncooperative vs. cooperative relax- 
ation. Like (i) and (ii), (iii) and (iv) are dichotomistic 
views of the character of the dominant relaxation mech- 
anisms. The "poles" correspond to whether the con- 
stituents of a system can be accurately treated as relax- 
ing independently of one another. If relaxation is unary 
and/or noncooperative, they can. If it is binary and/or 
cooperative, they cannot; instead, one must explicitly 
treat correlations between constituents. In systems with 
noncooperative relaxation, a tracer particle would have 
the same motion in a system where all other particles are 
frozen as in its unfrozen counterpart [1] . In systems with 
cooperative relaxation, the opposite is true. 

An example illustrating the unary vs. binary di- 
chotomy is as follows; suppose the stress u is a function of 
the intra-chain statistics (i?^(n)) = \ fi — fi+„|^, and per- 
haps the history of (i?^(n)), where n is chemical distance. 
In other words, suppose the stress is controlled by chain 
configuration but interchain correlations are either unim- 
portant or trivially integrable. Then the stress relaxation 
processes are unary. The classical theory of rubber elas- 
ticity, which assumes conformations of single strands at 
the scale n = Nc (the distance between crosslinks) con- 
trol stress, is a unary-relaxation theory. On the other 
hand, suppose local interchain orientation is important; 
relaxation processes will then be binary or higher order. 

A paradigmatic higher-order process in polymers is dis- 
entanglement, which occurs when one chain end deinter- 
sects another chain. This is a binary process because 
it involves two chains. While it is often useful to ap- 
proximate disentanglement as an "infinite" order process 
and "wrap" it into a unary mean-field theory (z. e., the 
tube theory of melt dynamics [3|), it is well known that 
in many cases such a description becomes inadequate Q . 
In principle, binary processes require a formal description 
that utilizes some two point correlation function F{fi, r2) 
where ri and r2 lie on different chains. 

Dichotomies (iii-iv) are are not new concepts, and 
have previously been applied in polymer melt rheology 
as well as many other fields. They are subject to 
ambiguities similar to those mentioned above for (i-ii); 
the same system may have unary and binary and / or 



cooperative and noncooperative relaxation processes oc- 
curring simultaneously on different length scales. Many 
processes are known to become increasingly cooperative 
as Tg is approached from above fl^ , and remain so below 
Tg, but predicting the degree of cooperativity is difficult. 
Another useful measure of cooperativity is the degree 
of coupling between different relaxation processes. This 
clearly strengthens as separation of the relevant EL&T 
scales decreases, but quantitative prediction is challeng- 
ing. Furthermore, in systems exhibiting nonlinear re- 
sponse, it is typically a priori unclear which couplings 
are most important in which regimes. Such issues remain 
controversial even for simple (e. g. metallic and colloidal) 
glassformers. However, since few theories of glassy poly- 
mer mechanics explicitly consider them, opportunities for 
progress abound. I now give an example illustrating how 
even qualitative treatment of dichotomies (iii-iv) can elu- 
cidate the nature of the crossover to dramatic hardening. 

Understanding the role played by entanglements in 
controlling the mechanical properties of polym er g lasses 
has proven extremely difficult. The canonical 111!] func- 
tional form of stress-strain curves at large strains, "Gaus- 
sian" strain hardening, is given by 

o- = (7o + G'flff(A), (1) 

where co is a flow stress, Gr is the strain hardening mod- 
ulus, A is the macroscopic stretch tensor, and 5(A) is 
alternately the (negative) derivative of the entropy den- 
sity of an affinely stretched entanglement network, or a 
Neohookean term. The traditional entropic treatment 
of strain hardening in glasses ^ qualitatively captures 
the shape of stress-strain curves, but has several flaws, 
extensively discussed in the literature (e. g. [H-[l3|)- 

Much progress has been made in the five years since 
Kramer challenged [l3| the polymer physics community 
to resolve this issue. Most physically, the mechanisms un- 
derlying strain hardening in glasses have been shown to 
be largely viscous and closely connected to plastic flow 
[13, Yet the underlying problem of distinguishing 

between models which assign different roles to entan- 
glements yet make similar predictions remains unsolved. 
Theories which assume entanglements are all-important, 
and theories which assume they play no role, both predict 
the Gaussian form. It is seemingly predicted by practi- 
cally any theory which predicts chains deform affinely on 
large scales, whether based on linearized entropic elastic- 
ity, Neohookean viscoplasticity [l^ or alteration of in- 
terchain ordering and suppression of density fluctuations 
0. Similarly, sub-Gaussian hardening, i. e. a sublinear 
in 5, is produced by subaffine large-scale deformation, 
whether this arises from relaxation of the entanglement 
network (as in melts or transient networks) or finite er 
in uncrosslinked glasses [12, 16]. Dramatic hardening (ct 
supralinear in g) can be produced either by entropic de- 
pletion of configurations for finite-n chain segments f3| or 
by the increased plastic deformation and bond stretching 
required to deform chains affinely while maintaining their 
connectivity. 
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FIG. 2: Crossover from unary-noncooperetative to binary-cooperative relaxation in polymer glasses arising from stretching of 
entanglement network. Systems are uniaxially compressed to a true strain e = —1.5.: A = exp(e), |(7(A)| = |A~^ — A^|, and 
— 1.5)1 — 4.43. Simulations are performed at low temperature T ~ Tg/35 to minimize thermal noise, and a high strain 
rate |e| = 10~^/tlj is purposefully chosen [17| . All units and protocols are desribed in the Appendix. Panels (a-b) contrast 
results from loosely and tightly entangled systems for (a) scaled stress a /{go + GRg{X)) and (b) energy dissipated (in units of 
Mo) per 'damaged' LJ bond. Statistical noise in panel (b) arises from finite system-size effects. Panels (c-d) show results for 
tightly entangled systems: variations with strain of (c) bond tension correlation fluctuations along chain backbones, and (d) 
the probability distribution for nonaffine jump sizes. 



Formulation of a robust microscopic theory predict- 
ing the entire range of mechanical response shown in 
Fig. [T] seems doubtful while these ambiguities remain. 
One approach to resolving them is to connect changes in 
macroscopic mechanical response with increasing strain 
to changes in relationships between structural features 
and relaxation mechanisms at different length scales. 
Figure [5] presents results for the crossover to dramatic 
strain hardening, which I will argue can best be under- 
stood as representing crossovers from unary to binary and 
from noncooperative to cooperative relaxation phenom- 
ena. Results are obtained via molecular dynamics simu- 
lations of a simple coarse-grained bead-spring model [l^ 
that captures the key physics of linear homopolymers. 
The simulation protocol is standard fllj and is described 
in the Appendix. Both flexible and semiflexible chains 
are studied to illustrate the behavior of "loosely" and 
"tightly" [20] entangled systems, respectively. 

Panel (a) compares the normalized stresses {a /{uq + 
Gfig^X))), where Gr is fit to the initial hardening regime 
(.5 < \g[X)\ < 1). The plateau for flexible chains in- 
dicates nearly Gaussian hardening at all strains. In 
contrast, semiflexible chains show dramatically supra- 
Gaussian hardening for \g{X)\ ^ 2. Panel (b) shows the 



energy dissipated per damaged bond Ud for the same 
systems. Ud = a'^/P, where cr'^ is the dissipative com- 
ponent of the stress and P is the rate of bond damage 
per unit strain (see the Appendix). Flexible chains again 
show a plateau, indicating Ud is constant when hardening 
is Gaussian, while semiflexible chains show an increase in 
Ud for (?(A) > 3. Differences between looselj^and tightly 
entangled systems are directly associated [14] with dif- 
fering degrees of increase in the energetic component of 
stress, = a — . 

Panel (c) shows the correlation in bond tensions T — 
dUpENE/di along chain backbones in tightly entangled 
systems. The correlations are roughly exponential. At 
intermediate n, {TiTi^n) exp{~2n/Ne) in the limit 
of large strains; the factor of 2 suggests a binary sup- 
pression of relaxation, wherein tension is concentrated 
at increasingly localized entanglement points (with 2 
chains/entanglement) and decorrelates between entan- 
glement points. In contrast, bond tensions in flexible 
systems (not shown) are much smaller and less corre- 
lated. This is of interest since chain tension relief by co- 
valent bond scission is a key mechanism leading to brittle 
fracture 

Panel (d) shows the probability distribution for dif- 
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ferential nonafHne jump sizes -D^„a in tightly entangled 
systems, defined by 



D2 

dna 



rk+i 



^k+l ^ 



(2) 



where fk is the position of a particle and is the macro- 
scopic stretch tensor at |e| = k6e. The tails of P{Ddna) 
become longer with increasing stress and strain [2l| . This 
effect is most pronounced for tightly entangled chains in 
the dramatic hardening regime (results for loosely entan- 
gled chains are similar to those found in Ref. [21| and 
are not shown). It arises because correlated bond ten- 
sion at the scale n ~ iVe increases the size of local plastic 
rearrangements, which in turn increases the energy they 
dissipate (panel b). 

All results in Fig. [5] are consistent with crossovers from 
unary to binary relaxation as chains stretch between en- 
tanglements, and from noncooperativity to cooperativ- 
ity between deformation at the level of the entanglement 
mesh and local plastic rearrangements at the monomer 
or Kuhn scale. For the strains considered here, these 
crossovers are present in tightly entangled but not loosely 
entangled systems; all coincide with increasing stretching 
of chains over chemical distances n ^ Ng. Note that the 
crossovers in Fig. [Ija-b) occur at different strains, while 
those in panels (c-d) are gradual. It seems likely that the 
unary-binary and noncooperative-cooperative crossovers 
are themselves coupled, though the strength of the cou- 
pling remains unclear. Understanding such behavior at 
a predictive level may be difficult, but presents an inter- 
esting challenge for the community. 

Single-chain-in-mean-field descriptions of strain hard- 
ening 0, @, [l^ probably cannot quantitatively treat 
the "entanglement-stretching" unary-binary crossover, 
even if interchain correlations are integrated into the 
mean field. Similarly, theories assuming a single relax- 
ation mechanism (e. g. segmental relxn.) cannot treat 
crossovers from noncooperative to cooperative relaxation 
which are driven by coupling to another relaxation mech- 
anism of different character. Presently, the microscopic 
theory which most satisfactorily captures the elastic, 
yield, softening, flow, and hardening regimes is due to 
Chen and Schweizer it is based on a strain- and 
thermal-history-dependent dynamical free energy for seg- 
mental rearrangements in the glassy state. While it cap- 
tures much of the physics of strain hardening (e. g. its 
coupling to plastic flow) by predicting how segmental re- 
laxation is suppressed due to changes in interchain or- 
dering in a macroscopic strain field, and even predicts 
variations in the character of activated segmental hops, 
it is questionable whether it can capture the crossovers 
shown in Fig. [2Kb-c), and it does not treat fracture. 

An interesting feature of the unary vs. binary di- 
chotomy is that unary relaxation processes "scale" lin- 
early with the density p of relaxing constituents, while 
binary processes scale quadratically with p. For exam- 
ple, the reason polymeric entanglement is traditionally 



regarded as a (nearly) binary process [23 is that en- 
tanglement density scales roughly as p"^, where p is the 
density of uncrossable chain contours. Studies of bidis- 
perse systems [1^ (with constituent densities (pi, P2) and 
relaxation processes which scale as P1P2) should there- 
fore be particularly useful in testing whether relaxation 
is unary or binary. Many such studies have been per- 
formed for melts well above Tg, in which the coupling 
between relaxation of short and long chains is now well 
understood , 24j , but similar studies in systems below Tg 
are in their infancy. New experimental techniques such 
as scanning near field optical microscopy [2^] should be 
particularly useful for these purposes. 

A specific example of how such studies can eluci- 
date causal relationships for polymer mechanics is as fol- 
lows. Experiments performed on (essentially) monodis- 
perse systems remain contradictory on such basic ques- 
tions as to whether or not Gr is proportional to entan- 
glement density pe [HI, [12: ^M '^ clear that 
the constant of proportionality differs for chemically dif- 
ferent polymers. In monodisperse (but not bidisperse) 
systems, both Gr and entanglement density pe scale ap- 
proximately as [H, |2^ . This "macro- micro" ambigu- 
ity led to much controversy over whether scales fun- 
damentally with entanglement density or Kuhn length. 
For systems exhibiting Gaussian hardening, an apparent 
resolution of these issues has been obtained via recent 
simulation and experimental studies of bidisperse sys- 
tems [ill, [23^ that suggest structure at the Kuhn level 
is the more fundamental controlling factor. However, the 
ambiguity remains unresolved in the dramatic hardening 
regime, and the crossovers described above suggest that 
developing a theoretical description that captures the en- 
tire spectrum of behavior depicted in Fig. [T] will require 
a great deal more work. 

In this effort, many of the most fruitful ideas may come 
from glass transition physics [To| . which I believe has 
been underutilized to date in theories of glassy polymer 
mechanics (at least those which treat large strains), ex- 
cept that of Ref. [6|. It may also be profitable to em- 
ploy concepts developed in recent studies of other "soft" 
systems, such as colloidal and granular materials. One 
promising possibility is to determine whether and when 
the properties of polyiners are more like attractive glasses 
or repulsive glasses [1, [l^l . Other potentially useful con- 
cepts include jamming, soft modes, inherent structures 
and energy landscapes [28| . Much also may be gained by 
working to bridge microscopic theories with increasingly 
sophisticated new constitutive models; some especially 
promising models incorporate melt-like relaxation mech- 
anisms [29j . I hope that this work will help spur further 
effort along these lines. 

The results presented in Fig. [5] are a new analysis of 
earlier simulations 1J| done in collaboration with Mark 
O. Robbins. Stimulating discussions with Kenneth S. 
Schweizer, Leon Govaert, Mya Warren, and Mark Rob- 
bins, and support from NSF Award Nos. DMR-1006537 
and DMR-0454947 are gratefully acknowledged. 
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I. APPENDIX: SIMULATION PROTOCOL 

All beads have mass m and interact via the trun- 
cated and shifted Lennard- Jones potential U^jir) = 
4wo[(a/r)^^-(a/r)'^-(a/rc)^^ + (a/rc)^], where = 1.5a 
is the cutoff radius and Ulj{t) = for r > re- The 
unit of time is tlj = y^ma^/uQ. Each polymer chain 
contains N beads. Covalent bonds are modeled using 
the finitely extensible nonlinear elastic (FENE) poten- 
tial U{r) = -(l/2)(fcEg)ln(l - {r/Ro f), with Rq = 1.5a 
and k = SOuq/o^ [l^, and have variable length with 
the equilibrium value £o = 0.96a. A bending potential 
Ubend(Q) = kbendi^ — cosO) , whcrc is the angle be- 
tween consecutive covalent bond vectors along a chain, 
imparts variable chain stiffness. Two values of k{,end are 
employed: flexible chains with kbend — have an entan- 
glement length Ne ~ 85, and semiflexible chains with 



kbend = 2mo havc Ne ~ 22 [l^]. The systems are well 
entangled; N = 500 for flexible and 350 for semiflexible 
chains. Periodic boundaries are applied in all three direc- 
tions, with cell dimensions L^. Ly, and along the x, 
y, and z directions. Well-equilibrated melts are rapidly 
quenched to fcsT = O.OImq ^ Tg/35. Uniaxial compres- 
sion is applied along the z-direction at constant strain 
rate e = Lz/L^, and and Ly are varied to maintain 
zero stress along the transverse directions. 

Simulation results were obtained using LAMMPS [s^. 
Bond damage (Fig.[51Jb)) corresponds to local plastic re- 
arrangements, identified with changes in intermonomer 
neighbor distances greater than 20% over a strain inter- 
val \Se\ = .025) il4j]. Nonaffine jumps (Fig. [l^d) and Eq. 
[5]) are defined using the same \Se\. The energy dissipated 
per damaged bond is given in units of mq, and bond ten- 
sions (Fig.mc)) have units uo/a. 
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